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Dynamical spin susceptibility is calculated for the generalized ferromagnetic 
Kondo model which describes itinerant eg electrons interacting with localized t2g 
electrons with antiferromagnetic coupling. The calculations done in the mean field 
approximation show that the spin-wave spectrum of the system in ferromagnetic 
state has two branches, acoustic and optic ones. Self-energy corrections to the spec- 
trum are considered and the acoustic spin-wave damping is evaluated. 



1 Introduction 



Manganites of the perovskite structure of the form Ri^xB^MnO^ (where R are triva- 
lent rare-earth and B are divalent alcahne ions, respectively) and related compounds 
that present the phenomenon of "colossal" magnetoresistance (CMR) have recently 
attracted much attention both from the basic point of view and due to their potential 
application [|T|, 0. The large magnetoresistance occurs close to the metal-insulator 
and the paramagnetic- ferromagnetic transitions where the interplay of transport, 
magnetic and structural properties is of the great importance (see Q). The key 
elements of the manganese oxides are Mn ions. In the parent compound LaMnO^ 
the electronic configuration of Mri^^ is {t^gCg). In this configuration due to a strong 
intra-orbital Hund's coupling t\g electrons go into tightly bound d^y, dyz, d^z core 
states and make up an electrically inert core Heisenberg spins S of magnitude 3/2. 
The configuration is very stable and remains localized over the entire range of 
doping. 



In the undoped case, with one e„ electron per Mn ion, two e„ orbitals, d^. 



and dsz^-r'^ types, are splitted due to the Jahn- Teller effect. At low temperature Cg 
electrons occupy d3x2_j,2 and d^y2_^2 ordered alternately in ah plane with their spins 
aligned to the core spin by interorbital Hund's coupling. Due to the Goodenough- 
Kanamori rules it results in the A-type antiferromagnet (AFM) ground state 
(AFM vector Q=(0, 0, 0.5)) with spin S = 2 for LaMnO^. Upon doping with holes 
by substituting La with Sr or any other divalent ions system becomes ferromagnetic 
(FM) and conducting. The hopping between two Mn sites is maximal when the core 
spins are parallel and minimum when they are antiparallel.That results in effective 
ferromagnetic exchange between the nearest neighbor core spins and thus leads to 
the FM metallic ground state of doped compounds. This behavior is qualitatively 
well described within the framework of the double exchange (DE) mechanism (see 
0,10, [0]). At higher hole concentration, x > 0.5, a charge ordering for holes is 
observed, and at x = 1 an insulating G-type AFM state (with Q = (0.5, 0.5, 0.5)) 
takes place for CaMnO^ compound. Therefore to describe the experimentally ob- 
tained phase diagram ( see, for example, pl) one should take into account both the 
Heisenberg type of AFM exchange between the core t2g electrons and the strong 
Hund coupling between and Cg electrons (see, e.g. |P, |l^). These competing 
interactions could be responsible for a coexistence of AFM and FM states observed 



recently in neutron scattering experiments in {LaQ,2hP'^o.7b)o.7CaQ,'iMnO^ [|TT| and 



in the bilayer manganite Lai,2Sri,^Mn20^ [T^. Also a crossover from an ideal 
isotropic FM spin-wave behaviour at low temperature to a diffusive spin propaga- 
tion observed in LaQ jCaQ^^MnO^ [|T^ could be explained if one takes into account 
both the localized t2g spin {S = 3/2) and the itinerant Cg spin (a = 1/2) subsystems. 

In the present paper we study the spin dynamics in manganites within the gen- 
eralized ferromegnetic Kondo model (FKM) allowing for the AFM exchange inter- 
action between t2g spins. Unlike to the DE model (see, [|14|), where Jn/t — oo 
is considered and the system is treated as perfectly spin polarized with S' = 2, in 
our work both the fluctuation of the localized and itinerant spins are taking into 
account. However, we ignore in the present calculations a possible orbital ordering 
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which is very important in explaining different types of AFM ordering in the insulat- 
ing phases ||10|, |15[ but plays less essential role in the FM state considered here. To 
take into account strong Coulomb interaction between Cg electrons which excludes 
the double occupancy of Cg electrons for a lattice site we employ the Hubbard op- 
erator technique. The spectrum of spin waves in the FM state is calculated by 
employing equations of motion for the matrix Green function (GF) for the localized 
and itinerant spins. In the next Section the model and general formalism for the GF 
are presented. The spin-wave spectrum in a generalized mean field approximation 
(MFA) is calculated in Sec. 3 and self-energy corrections and spin-wave damping 
are evaluated in Sec. 4. 



2 The model 

We consider an effective Hamiltonian of the generalized FM Kondo model which 
can be written in the following form [Q: 

^ = - E U.XfX^'^ - g 5: S,a„ + 1 E J,,S,S, . (1) 

i,j,a i i,j 

The first term of Eq. describes an electron hopping between Mn-ions where X°"° 
is the creation operator of an electron with spin a in one of the eg orbitals. Here we 
neglect orbital degeneracy of eg electrons and introduce orbital independent hopping 
parameter tij with tij = t for the nearest neighbors. The second term describes the 
ferromagnetic Hund coupling {J^ > 0) between eg and spins where S, refers to 
the localized Mn core spin 5 = 3/2. The third term describes the antiferromagnetic 
coupling of localized spins between the nearest neighbor sites. In real materials the 
coupling of core spins is not the same in different directions and should be described 
in the matrix form, but for simplicity we are analyzing the isotropic case [Jij = J). 
We exclude the doubly occupied eg state from the effective Hamiltonian by using 
the Hubbard operator representation because the electron-electron interaction has 
the largest energy scale (intra-atomic Coulomb interaction in the eg orbitals) and 
can be estimated as 7 — 8 eV while Jh ^ ^ eV. Due to large Hund energy we 
neglect superexchange interaction between eg electrons of the order of t'^/U ||10| . 
The conduction bandwidth is smaller than the Hund coupling energy and from 
density-functional studies can be estimated as t ~ 0.15 eV |T6| . 

The HQ's in Eq. (|l]) are defined as Xf^ = \i,a){i,P\ for three possible states at 
the lattice site i: \i,a) = \i,0), \i,cr) for an empty site and for a singly occupied 
site with spin a = (|, |) = (+, — ) . The completeness relation for the HQ's reads as 

^r+E^r = i- (2) 

cr 

For itinerant electrons the spin and density operators in Eq. (|T|) are expressed by 
HQ's as 

at = Xl\ a- = Xl\ at = ^(X/^ - X^), n, = X}^ + X^. (3) 
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The HO's obey the following commutation relations 



aS 



5ij ((5/3^Xf ± SsaXi 



(4) 



In Eq.(^) the upper sign stands for the case when both HO's are Fermi-like ones 
(as, e. g., X^'^). The spin and density operators (H) are Bose-like and for them the 
lower sign in Eq.(|^) should be taken. 

It is assumed that the core spin operators Sf obey the standard commutation 
relations, e.g.. 



(5) 



To treat the fluctuations of localized and itinerant spins at the same level of 
approximation we introduce the dynamic spin susceptibility (DSS) of the system in 
the matrix form 

f Xn Xi2 ^ 

\ X21 X22 , 



xiq,^) 



(6) 



where 



Here 



A. 



5^ 



4 + 



/■oo 1 



(7) 



denotes the Fourier transformed two-time retarded commutator Green function (GF) 



Tq , [T9|| . The diagonal elements xiiiQ,^) and X22{(1,^) stands for the itinerant and 



core spin GF, respectively, while the nondiagonal elements XuiQ,^) and X2i{(l,^) 
define the crosscorrelations between the two spin subsystems. The GF (P) obeys the 
following equation of motion 



u{{^A,\A^))^ = + {{^A,\-^Ap)^ . 



These equations (j^) could be easily combined in a more convenient form of the 
equation of motion |]T^: 

u{{A,\A^))^ = {[A„A^]) 



{[A,,Aj]) 



■ {{AMt)l 



(9) 



where the current is defined as tA = idA/dt = [A,H] and in the matrix form can 
be given by the following expression: 



tA„ 



(10) 



and 



0)r 



- - (11) 

is the irreducible part of the higher order GF. 
We can rewrite in the Dyson form 

Xgico) = -n,- n(g, to)]-' ■ I, (12) 

where Tq is the unity matrix and 

where {a') = {af) and {S') = (5f) . 

The matrix Qq = Qq describes the mean field (MF) energy spectrum and 
Il{q,uj) = Il{q,uj) /-^ is the self-energy matrix. They are given by 



with 



^^r = E^^K^f^r - xrx^') - ^{sta! - Stat) , (16) 



2S 



^St = -^{S[at - Stat) - E JaiStSt - S^Sf) . (17) 



3 Mean field approximation 

Let us now examine the spectrum and DSS in mean field approximation (MFA). 
The spin-wave dispersion is determined by the following equation 



det [uotq - figj = . (18) 
From ( [I^ ) we obtain for the matrix elements of Vlq 

~ _( [d + a{l-^q)]/2{a^) ~d/2{S^) \ 

' [ -d/2{a^) [d-b{l-^q)]/2{S^) ) ' ^'^^ 

where we are using the following notation: 

d=^{2{afSn + {atSr)), (20) 
a = zt{nl + n{), b = zJNi, (21) 



5 



with tq = ztjq, 7g = {2/z) (cos Qx + COS Qy + COS Qz) , whcie 2; = 6 for the simple three- 
dimensional cubic lattice with nearest-neighbor hopping t. In ( ^1]) the nearest 
neighbor particle-hole and spin correlation functions are defined as follows 



N 



Til. 



(22) 



The equation (|TB|) has two solutions describing two branches of spin wave excitations: 

(23) 



Kl(2) = t 

^ 2 
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For the model calculation we can expand this equation at g and for the finite 
value of d we obtain 



Ef) ~ A + D^q^ 



(24) 



where E^^ corresponds to the gapless (acoustic) spin-wave excitation with the stiff- 
ness Di given by 

a — h 

' ^ i2(MTM) ' 

and Ejj^^ describes the optic mode of the spin fluctuations with the gap A and the 
effective stiffness D2 determined by the following expressions: 



A = d 



Do 



(26) 



The ferromagnetic acoustic spin-wave becomes unstable when the stiffness Di ^ 
or a — 6 = in Eq. (p5l). It may happen for small concentration of itinerant electrons, 
n < ric — 2SJ/t ~ 0.3. The self-energy corrections considered below (see Eq. (^)) 
even increase the critical value ric- 

The spectrum of spin fluctuations in MFA are given by the spectral functions 



(27) 



for the spin susceptibility 



In 



22 



E, 



(2) 



-6{u;-E(''>) + 



1 a 



E, 



(2) 



f28) 
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B^/ig,u) = 2{S^) ;,,, 6iu - Ejp) + ^ J,, 6{u - )) , (29) 

\rjq — rbq tjq — tjq I 

Br/{qM = B^/{q,u) = ^ (5(^ - i?«) - 5(^ - i^f )) . (30) 

-C/g — Eq 

The spectral functions (|28|) - ( pOD obey the following sum rules: 

/ + 00 
-oo 

/-l-oo 
u;c?u;S^J^(g,u;) = fif = (M„ , (31) 

where the matrices 1^/3 and are given by the Eqs. ([T3|) , ([T^) . 

4 Self-energy corrections 

The next step is to consider the self-energy corrections to the MF spectrum. Taking 
into account the self-energy corrections the equation for the spectrum transforms 
into the following form: 

det (euro - fig - fl{q, cu)) = 0. (32) 

First we compute the self-energy matrix elements by using mode-coupling approxi- 
mation in terms of the dressed particle-hole and spin fluctuations (see, e. g., Gotze 
et al., [|l7l)- This scheme is essentially equivalent to the self-consistent Born ap- 



proximation in which the vertex corrections are neglected. The proposed scheme is 
defined by the following decoupling of the time-dependent correlation functions: 

{at{t)S;^{t)apr) ^ {at{t)a^^){S;MSt) , (34) 

{St{t)s;^(t)s^^st) ^ {St{t)s^^){s;^{t)st) . (35) 

The self-energy matrix elements are obtained by using the above defined decou- 



pling scheme (|33|) and (^) with the spectral representation for the GF. The diagonal 



elements involve two contributions: 

nii(22) (g, w) = ^u(22) ^) + nS?(22) ^) • (36) 

The first one describes fluctuations of the internal degrees of freedom of the given 
spin subsystem. While the second one stems from the Hund's term and describes 
the coupling between itinerant and core spins. 

For the itinerant spins the first term in Eq. ( |5BD is due to decay of spin fluctua- 
tions into particle-hole pair excitations and reads as 

1 r+co r+oo n{uji - uj') - n{uji) 



nil (?,^) = T7 / / duj'dui 

i V J — oo J — oo 



UJ — uj' + IS 
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X E tlA'^ik - g, ^1 - uj')A^k, 00,) (37) 

k,a- 

where tkq = zt{jk — Ik-q), n{uj) = (e^'^ + l)^"*^, and A'^{k,u!) is the single-electron 
spectral function. By using the MF approximation for that, Eq. (A. 7), we can 
integrate over the frequencies in Eq. ( ^7|) and obtain the following estimation: 



^niq,u;) = 1 (1 - - n'^) ^^"^'^^ ""^''^ . (38) 

It has the standard form for a one-loop particle-hole contribution to the self-energy 
(see, e.g. 0). 

The second terms in Eq. (pG]) are the same for both subsystem and coincides with 
the nondiagonal elements of the self-energy matrix due to the Hund coupling 

nS?(22) ^) = -ni2(2i) (g, ^) = (g, w) (39) 

with 

\2S J Ntx'^ J-oo J -00 uj — oj' + ie 

X E [IniX22(^, uji)luixii{.k -q,uj' - uJi) + Imx22(A;, uji)lmxliik - q,uj' - Ui)] (40) 

k 

where N{u) = (e^"^ — 1)^^, Xu X22 denotes the longitudinal susceptibility of the 
itinerant and core spins, respectively. 

Let us consider now the remaining 1122^ term which describes the fluctuations 
in the core spin subsystem. This contribution is due to the Heisenberg exchange 
between the localized spins and is given by 



n22(g,-) = ^y_ y_ 



1 + N{uj' - LJi) + N{ui) 



X 



00 J -00 uj — uj' + te 

E Jkq^^X22{k - Q, UJl)lmX22{k, Uj' - UJl) , (41) 



where Jkq = zJ{-fk - Ik-q)- 

In order to evaluate the longitudinal susceptibility in Eqs. (|3)' (0) ^^r both 
subsystems we will use for them the simplest one- loop approximation (see, e.g. @). 
In this approximation the imaginary part of Xiiil)^) is given as the convolution of 
the single-electron GFs 

\ \ r+co 

ImXii(g, ^) = 7T7 / duj'[n{uj' - uj) - n{^')] 

xJ2A''{k,uj')A^k-q,uj' -00). (42) 

k,(T 

The imaginary part of the core spin susceptibility X22 can be expressed in the linear 
spin-wave approximation as 

--Imx^2(g,^) = ^^^^ y_ du;\N{u;'-u)~N{u;')) 
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J2 Imx22(fc -q,^' - uj)lmx22{k, uj') , (43) 

k 

which follows directly from the Holstein-Primakoff representation. 

To study the spin wave spectrum including self-energy corrections let us consider 
the static limit for g ^ 0. For the self-energy matrix we can write 

where 

A = - M,„ B = - 1„„ 5l(|l£) , 4 = -n^»(0.0). (45) 

Here the coefficient A, B, and di are positive since n(g, 0) < in the second 
order of the perturbation theory. As it follows from Eqs. (|^) and (^), the self- 
energy corrections coming from the Hund's coupling does not renormalize the spin 
stiffness and gives input only into the gap. Hence the spin-wave spectrum in the 
longwavelegth limit can be written as 

cjf ^ A + D^q^ (46) 

where the renormalized spin stiffness and the gap are given by 

a-b-A-B 



12{{S^) + {a^)) 
~ ^ {a-A){S^)'-ib + B){ar 

Now we consider the damping for the acoustic spin wave mode given by the 
imaginary parts of the self-energy matrix. For the damping induced by particle-hole 
excitations we get from Eq. (^8]) in the MF approximation for the single-electron 
GFs: 

rSV(g,cu) = -ImHff(g,a; + 2£) 
= ^ E4(l - - ^nHel-q) - nielMu + el_^ - 4) (48) 

fc,(T 

The contribution, due to the finite /c— independent gap in single-electron spectrum 
in the ferromagnetic state disappears in the low frequency limit, T^ii{q^uj) = for 
00 <h (see Eq. (A. 14)). 

The damping due to the antiferomagnetic exchange interaction given by the 
imaginary part of the self-energy Il22{q,^^), Eq. (0), gives a small contribution 
proportional to in the longwavelength limit and can be disregarded due to small 
antiferromagnetic exchange interaction J. 
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The largest contribution to the damping of spin waves is given by the imaginary 
part of the self-energy due to the Hund coupling, Eq. (^): 

/ 7„\ 2 1 r+°° 

Tniq, uj) = -Imn^^(g, lu + te) = i^^j — J^^ du.il + N{uj - Ui) + N{uj,)) 
X [ImX22(A; - - uji)lmxuik,uji) + Imxii(A; -q,uj- u;i)Imx22(^, ^i)] • (49) 

k 

Here for the imaginary parts of the spin susceptibilities Xiii^^^) X22{k,u!) we 
will use their MP values in Eqs. (^), ( PU[ ) taking into account only the acoustic E^^'' 
mode: 

- -lmxr/{q,u;) ^ 2(a^)-|^-^<5(^ - Ejp) = Af6{u - E,) , (50) 

Eq — Eq 

1 9}^ - K(i) 

- -Imxr(g,a;) ^ 2(g-) J ^ 5{uj - i?«) = Af 5(c. - Eq) . (51) 

^ Eq — Eq 

For the longitudinal spin susceptibilities in (^) we will use Eqs. (0), (^). 
After integration over uji we get the following result 

Vn{q.u) = (e-/- - 1) (g)' ^ E/f-.(l - n.f 



xNiE,^q)niel)[l - n(4^_,)]5(^ - E,.q + e^^, - el) 

^(p'^/T \ f:lli\ ^ V All A 22* 22 

) \ OQ AQ2AT2 2^ ^^kW^ki^^ki-k 



\2SJ Asm\,^^ 

X NiE,.q)NiEkJ[l + NiE,,,,)]6iuj - Ef^l^ + E^^l, - E,J . (52) 

It describes a spin wave damping due to its decay into an electron-hole pair and 
another spin wave (the first term) and a three spin-wave scattering process (the 
second term). The latter has a standard form for three magnon scattering (see, 
(31.2.20) in |^). At low energy regime {u — >■ 0) and at the longwavelenght limit 



the requirements for conservation of energy and momentum allow only small wave 
vectors and thus only small energies. Hence we can consider A^^ and A|^ as k- 
independent. In the limit A; — we obtain 

. a;- ^ 2W')^^^ , Af . AS' ^ 2iS')j^^ . (53, 



In this approximation from the Eq. ( |52D follows that at low energies, (cu <^T), the 
damping has a linear tu-dependence, T'^"{q,uj) ~ u, and does not depend explicitly 
on the wave vector q. In the low temperature region, (a; <^ T -C cuq) where c^o 
is the maximal acoustic spin-wave frequency, estimations for rH{q,uj)/u! show 
that the first contribution in Eq. (^) is proportional to {T / Uq) {Jjj / N {eF)vpko) 
where N{eF) is the density of state at the Fermi level and vp is the Fermi velocity, 
ko ~ 27r/a and a is a lattice constant. The second term gives to the damping the 
contribution proportional to {T /ujq)[Jh/<jJo)'^ ■ To give more accurate estimations 
for the spin-wave damping numerical studies should be performed which will be 
considered elsewhere. 
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5 Conclusions 



In the present paper we have calculated dynamical spin susceptibility for the gen- 
eralized ferromagnetic Kondo model (1) by taking into account explicitly both the 
strong Hund interaction for the itinerant Cg and localized t2g electrons and AFM 
interaction between the t2g electrons. We consider the ferromagnetic phase and 
therefore neglect a possible orbital ordering of the Cg electrons. Strong electron cor- 
relations between Cg electrons are treated within the Hubbard operator technique 
which is important in calculation of the single-electron GF for the itinerant electrons. 

We have proved that even in the MFA described by the frequency matrix, 
Eq. (19), we get the , acoustic spin-wave excitations, Eq. (24), due to coupling 
of the two modes with gaps for itinerant and localized electrons. The gapless mode 
should appear in the model (1) with rotation symmetry for spin system. A gapless 
mode in the limit Jh —>■ oo, considered in Ref. |]l3, was obtained only by taking into 



account self-energy corrections. In our case the self energy corrections calculated in 
the self-consistent Born approximation, Eqs. (32)-(34), resulted in additional renor- 
malization of the stiffness of the acoustic ferromagnetic spin waves. The imaginary 
parts of the self-energy gives the damping of spin waves. We have evaluated the most 
important contribution due to Hund coupling in the second order, Eq. (49), which 
can be described as a three magnon scattering. The damping of acoustic spin waves 
is proportional to the frequency for u <^T, Eq. (51), and should be small for small 
wave vectors. To give numerical estimations one should solve self-consistently the 
system of equations for the matrix spin susceptibility, Eq. (6), and the self-energy 
functions, Eq. (37). Also the spectrum of single electron excitations, given in the 
Appendix, should be evaluated to consider the itinerant electron contributions to 
the spin waves, Eq. (36). These problems will be considered elsewhere. 
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Appendix 

In this Appendix we evaluate the single-electron Green function defined as 

G'^{k,u) = {{X^'^\Xf))^ (Al) 

with the corresponding spectral function Al{uj) = — (l/7r)Im G"{k^uj + le). In the 
site representation the equation of motion for G^{k,uj) reads as 

uj{{Xr\Xf))^ = ({xf^Xf }) + ((zX^lXf)). . (A2) 

Using the commutation relations, Eq.(^), we obtain 

txr = h] = -Y: tu [(xf + xDx^ + xrx^^' 
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■Jh 

'as 



The next step is to define the irreducible part of the current operator iXf'^ by 

(AA) 



{A.3) 



I 

The definition gives for the frequency matrix 

where by using the completeness relation, Eq.(|^), we write {^X^'^ , Xj'^ 
rf) with = (n°"). By using equation of motion (A. 3) we get 

({zXf^Xf }) = -t,,((l - nf)(l - nf)) - ^{XrXf) 



{A.h) 



(A.6) 



In the present paper we will not include the self-energy correction coming from Z^" 
term (A. 4) and treat the single-electron GF within the linear, MF type approxima- 
tion. That results in the following form of the single-electron GF 

1 _ r)0 

G"^(A;,cu) 



By introducing the nearest neighbor charge-spin correlation function 

1 



(A.7) 



N 

= {xrx!::) + (xf^xf/j 

we can rewrite the frequency matrix in the form 
where 



t,,[(l-2n'^) + X-]/(l 



n 



is the one particle spectrum in the linear approximation, and 



e = - 



Jh_ 

AS 



a{S^) 



{SfXf) 
1 - ra<^ 



+ 



ztnl 



(A8) 

(A.9) 
(AlO) 

(All) 



is the spin- dependent energy shift of the spectrum. In the momentum space the 
spectrum is 

4 = E e-^'^'^er, = - <//7., (^-12) 

where 

tiff = t[(l - 2nn + X- ]/(l - (A.13) 

is an effective bandwidth that is narrowed and spin dependent due to the spin and 
charge correlations. In the MFA in the ferromagnetic state we have a spin gap in 
the single-electron spectrum: 

el-elc:^ a{JH/2S){S') = ah. (AAA) 
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